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ABSTRACT 

We report observations at millimeter and submillimeter wavelengths of comet 
C/2009 PI (Garradd) from 2011 December 28 to 2012 April 24, using the Arizona 
Radio Observatory submillimeter telescope (SMT) and the James Clerk Maxwell 
Telescope (JCMT). Garradd is a dynamically young long-period comet from the 
Oort Cloud, with a periodicity of 127,000 years, that reached perihelion on 2011 
December 23 (at Rh = 1.55 AU and A = 20.1 AU) and made its closest approach 
to the Earth on 2012 March 05 (at R h = 1.84 AU and A = 1.26 AU). We 
obtained gas production rates, and molecular abundances relative to water for 
HCN, ortho- H 2 CO, CS, CO and CH 3 OH. A rotational temperature, T rot ~ 50 
K, was determined by observing multiple methanol lines with the JCMT. By 
averaging the abundance ratio relative to water from the SMT and the JCMT 
we derive: CO: 7.03 ± 1.84 %, HCN: 0.04 ± 0.01 %, o-H 2 CO: 0.14 ± 0.03 % as a 
parent molecule (and 0.28 ± 0.06 % as an extended source), CS: 0.03 ± 0.01 % and 
CH 3 OH: 3.1llg;5® %. We concluded that Garradd is normal in CH 3 OH, depleted 
in HCN, o-H 2 CO and CS and slightly enriched in CO with respect to typically 
observed cometary mixing ratios. We also studied the temporal evolution of HCN 
and CO and find that the production of HCN has a trend similar to water (but 
with short-term variation), with a decrease after perihelion, while that of CO 
shows contrary behavior: remaining constant or increasing after perihelion. 

Subject headings: Astrobiology - Solar System: general - Comets: individual 
(C/2009 PI (Garradd)) - techniques: spectroscopic - submm lines: Solar System 
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1. Introduction 


Comets are among the most pristine objects in our Solar System. Their chemical com- 
position provides important clues to the processes that occurred during the formation and 
early evolution of the Solar System. Cometary nuclei today reside in (at least) two distinct 
reservoirs, the Oort Cloud (OC) and the Edgeworth-Kuiper Belt (EKB) (divided into the 
classical KB, the scattered disk, and the detached or extended disk populations; Morbidelli 
et al. (2008)). Once injected into the inner planetary system, comets are classified dynami- 
cally as nearly isotropic [long-period comet (LPC) or ffalley-type comet (HTC)] or ecliptic 
[Centaur-type, Encke-type, or Jupiter-family comet (JFC)]. The Tisserand parameter of the 
current orbit of a given comet identifies the storage reservoir from whence it came; eclip- 
tic comets come from the scattered KB reservoir, whereas the storage reservoir for nearly 
isotropic comets (NICs) is the OC (Levison et al. 1996; Horner et al. 2003). Past observations 
have shown that comets appear to contain a mixture of products from both interstellar and 
nebular chemistries and could also have been important for initiating prebiotic chemistry 
on the early Earth (Ehrenfreund & Charnley 2000). Although there are some differences, 
the volatile composition of cometary ices is generally similar to the inventory of molecules 
detected in the ices and gas of dense molecular clouds. Given the gradient in physical con- 
ditions expected across the proto-Solar nebula, chemical diversity in the comet population 
is to be expected, as has been inferred for both JFCs (Crovisier et al. 2009; Gicquel et al. 
2014) and OCs comets (DiSanti & Alumina 2008). 

Comet C/2009 PI (Garradd), hereafter Garradd, was discovered on 2009 August 13 
at a heliocentric distance Rh = 8.7 AU (McNaught & Garradd 2009). The comet reached 
perihelion on 2011 December 23 at Rh = 1.55 AU and perigee at A = 1.27 AU on 2012 
March 05. Observations were made throughout this period by many groups using both 
ground-based and space-based telescopes (e.g., Paganini et al. (2012), Bodewits et al. (2014) 
and Villanueva et al. (2012)). Based on the inclination to the ecliptic, Tisserand parameter, 
and original semi-major axis, comet Garradd can be classified as a probable dynamically 
young, long period comet from the Oort Cloud, whereas a classification based on its activity 
is less definitive, due to the diversity in composition both among and within the different 
comet populations. 

New dynamically young comets from the Oort Cloud (on their Erst trip back to the 
inner Solar System) may reveal primordial composition while the composition of comets in 
short orbits may be altered with each passage. Also, probably due to a different composition 
or the effect of the Solar heating (Oort & Schmidt 1951; A’Hearn et al. 1995; Cordiner et al. 
2014; Combi et al. 2014), dynamically new comets (like ISON) have shown different behavior. 
However, they tend to be more active during their early approach to the Sun and develop 



3 


more normal activity later in the apparition (Meech & Svoren 2004; Meech et al. 2009). 

Garradd’s pre-perihelion peak in water production between 100 - 50 days before peri- 
helion (Combi et al. 2013; Bodewits et al. 2014) and monotonically increase of CO activity 
post-perihelion (Feaga et al. 2014) is in contrast with the behaviour of the well studied 
Oort-Cloud comet C/1995 01 (Hale-Bopp), which showed an increase of all gas production 
rates pre-perihelion and a decrease afterwards (Biver et al. 2002). Garradd has also been 
defined as CO-rich, having a abundance ratio of CO relative to water between 9.12 ± 0.80 
and 13.45 ± 1.45 % (Paganini et al. 2012; Villanueva et al. 2012; DiSanti et al. 2014), which 
is consistent with other active Oort Cloud comets like C/1995 01 (Hale-Bopp) and C/1996 
B2 (Hyakutake). 

Intriguing early observations of comet Garradd, presented an opportunity to expand 
the still small database of observations of comets, from which eventually will emerge an 
understanding of their formation, evolution, and perhaps taxonomy. 

We conducted observations at submillimeter wavelengths of comet Garradd from 2011 
December 28 to 2012 April 24. Here we report detections of HCN, o-H 2 CO, CS, CO and 
CH3OH using the Arizona Radio Observatory submillimeter telescope (SMT) and the James 
Clerk Maxwell Telescope (JCMT). We present the observational results and determine phys- 
ical parameters such as rotational temperature, column densities, production rates, and 
abundance ratios. 


2. Observations 

All the observations presented here were obtained using the JPL/Horizons (Jet Propul- 
sion Laboratory) ephemeris number 82. 


2.1. Arizona Radio Observatory Submillimeter Telescope (SMT) 

Observations of HCN:J=3-2 (265.8864 GHz), CO:J=2-l (230.5380 GHz), CS:J=5-4 
(244.9356 GHz), o-H 2 CO:JK a ,K c =3i l2 -2i,i (225.6978 GHz), and CH 3 OH:J=5-4, toward comet 
C/2009 PI (Garradd) were conducted between 2011 December 28 and 2012 March 06 using 
the 10m Arizona Radio Observatory Submillimeter Telescope (SMT) on Mount Graham, 
Arizona. These data were obtained with a dual-polarization ALMA Band 6 receiver system, 
employing sideband-separating mixers with an image rejection of typically 14-25 dB. The 
backends employed were a 2048 channel 1 MHz and 250 kHz filter banks used in parallel 
mode. The temperature scale at the SMT is Tjj; radiation temperature is then defined as 
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T r = T* A /rj bl where r)b is the main beam efficiency. Data were taken in position-switching 
mode with an off position 30’ west in azimuth. The pointing accuracy is estimated at 1” 
rms. Focus and positional accuracy were checked periodically on nearby planets and masers. 


2.2. James Clerk Maxwell Telescope (JCMT) 

Observations of C/2009 PI (Garradd) were made from the JCMT, located at the 4000m 
level on Mauna Kea, Hawaii, on multiple occasions from 2011 December 30 to 2012 April 24, 
using the HARP heterodyne array in single-receptor mode. For this work, HARP was tuned 
for the transitions of HCN:J=4-3 (354.5055GHz) and CO:J=3-2 (345.7986GHz), employing 
the ACSIS digital correlation spectrometer configured at its highest frequency resolution (31 
kHz), and subsequently smoothed prior to presentation here. 

At these frequencies, the JCMT has a Gaussian beam of size « 15” (full width at half 
power). Pointing of the telescope was verified approximately every hour and was achieved by 
doing spectral-line ’five-points’ on nearby line sources. The pointing accuracy is estimated 
at 2” rms in each of two orthogonal coordinates (azimuth and elevation). Telescope focus is 
similarly maintained throughout the night by measures of bright line-sources. Calibration 
of the (T* a ) brightness scale is achieved by making measures at standard frequencies (e.g. 
CO:J=3-2, 345GHz) of astronomical sources used as reference calibrators. As for the SMT, 
the radiation temperature is defined as Tr = TX/v b - 

The opacity of the sky above JCMT was measured by a water vapor meter (WVM) 
mounted so as to measure along the telescope line-of-sight. Opacity is expressed as if mea- 
sured at the zenith at 225GHz. The opacity (in nepers) on UT 2011 December 30, 2012 
January 20 and 22, February 29, March 26, and April 24, was on average, 0.05, 0.10, 0.06, 
0.14, 0.18, 0.07, respectively. Spectroscopic data were reduced using STARLINK software 
smurf /makecube and analyzed using kappa/splat. 


3. Results 

Circumstances of the observations - UT date, observing frequency (u), beam size (© 5 ), 
diameter of the projected beam size on the comet (D), main beam efficiency (rjb), heliocentric 
distance (Rh) and comet distance at the times of measurements (A) - are listed in Tablet. 
The integrated line intensity (f T* A d v ) derived from a Gaussian fit and the associated la 
uncertainty are also given in Table. 1. 

Representative spectra from both facilities are shown in Figure 1, plotted in a come- 
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Table 1. Observations HCN, CO, CS, o-H 2 CO and CH 3 OH toward comet C/2009 PI 

(Garradd) from SMT and JCMT. 
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Reference 
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HCN 

SMT 

3-2 

2011 

Dec 28.94 

a 
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28.4 

40568 

0.74 

0.32 

± 

0.02 

1.55 

1.97 




2011 

Dec 29.33 

b 



40465 


0.26 

± 

0.02 

1.55 

1.96 




2012 

Jan 18.52 

c 



35008 


0.30 

± 

0.02 

1.59 

1.74 




2012 

Jan 19.53 

d 



35481 


0.35 

± 

0.02 

1.59 

1.72 




2012 

Feb 18.46 

e 



27739 


0.17 

± 

0.01 

1.73 

1.35 




2012 

Feb 19.36 

f 
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0.31 

± 

0.02 

1.74 

1.34 
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Feb 23.62 

g 



26894 


0.26 
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0.01 

1.76 

1.31 
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h 
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0.28 
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0.02 

1.84 

1.27 

CO 

SMT 

2-1 

2012 

Jan 18.63 
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0.04 

± 
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0.01 

1.73 
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1 
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0.02 

2.25 

2.02 



6 


tocentric velocity frame, and labelled by a letter corresponding to the entries in Table 
1. Figure 1 shows the CO:J=3-2 and HCN:J=4-3 transitions observed from comet Gar- 
radd on 2012 March 26 from JCMT, and the CO:J=2-l, HCN:J=3-2, CS:J=5-4 and o- 
H 2 CO: Jx a ,A' c =3i,2 — 2i,i transitions observed on 2012 March 06 from the SMT. Gas pro- 
duction rates from these transitions are reported in Table 2. Figure 2 shows the spectra 
of CH 3 OH:J=5-4 taken on 2012 February 19, and March 26 with the SMT. Derivation of 
the rotational temperature of CH 3 OH can be found in section 4.1, the computation of the 
column densities, production rates and abundances ratios is in section 4.2, and a detailed 
discussion about the temporal evolution of water, CO and HCN is in section 4.3. All data 
are available in the online journal as supplemental information. 



Fig. 1.— Detection of CO (J=3-2), CO (J=2-l), HCN (J=4-3), and HCN (J=3-2), o-H 2 CO 
{JK a ,K c — 3i, 2 — 2i,i), and CS (J=5-4) taken on 2012 March 06 (k, h, n and 1) and March 26 
(u and w) with the JCMT and the SMT. See Table 1 for dates, parameters of observations, 
and reference letters on the right-hand side of each panel. Additional spectra for the 
observing campaign for CO, HCN and o-H 2 CO are available in the online journal. 
Note that the plots n, h, w, k and u are offset from the baseline at zero. 


4. Analysis 

4.1. Rotational temperature 

We have applied the rotational diagram technique (Schloerb et al. 1983; Blake et al. 1986; 
Bockelee-Morvan et ah 1994; Drahus et ah 2012) to our SMT methanol data to determine 
the coma temperature and the total number of methanol molecules observed in the beam. 
We assume local-thermodynamic-equilibrium (LTE), thus a linear ln(F)=f(E u ) provides the 
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temperature and the column density: 


In 


Nu\ 
9u ) 



E, 


up 


kX, 


rot 


(i) 


with : 


N u 


8nku 2 f T^Avi /2 
hc 3 A u i‘r/ b 


( 2 ) 


where N u [cm -2 ] is column density upper energy level, k is the Boltzmann constant, ( rot is the 
partition function, h is the Planck constant, c is the speed of the light, A u i [s -1 ] is the Einstein 
coefficient, g up is the statistical weight, T rot [K] is the rotational excitation temperature, E up 
[cm -1 ] is the upper state energy, and N to t [cm -2 ] is the total number of molecules observed 
in the beam. We observed methanol lines on 2012 February 19 ( Rh =1.74 AU) and March 
06 (Rh = 1.84 AU) with the SMT, see Figure 2. The intensities of the Eve and six observed 
methanol lines (Figure 2), allowed us to determine a poorly constrain rotational temperature 
Trot = 46lfg K (2012 Feb. 19) and T rot = 46^ K (2012 March 06), or T rot & 50 K. Our value 
is in good agreement with Biver et al. (2012), who derived an average of the gas temperature 
in the coma of ~50 K from observations of 70 lines of methanol obtained with the IRAM-30m 
telescope on 2011 October and 2012 February. Using JCMT data taken on 2012 July 28-31 
and January 06-08, Yang & Drahus (2012) showed that T rot varied with R h , going from T rot 
ps 30 K at Rh = 2.49 AU, to T rot ~ 40 K at Rh = 1.56 AU. From the HON lines observed 
on 2011 September 8 and 9 using the high-resolution infrared spectrometers (NIRSPEC at 
Keck II and CHSELL at IRTF), Villanueva et al. (2012) determined T rot = 40 ± 7 K (Rh 
= 2.1 AU). Also from infrared data (CRIRES at ESO’s Very Large Telescope), Paganini 
et al. (2012) concluded that a T rot ps 50 K is consistent with the average temperature from 
the volatiles of H 2 0, CO, C 2 H 6 , HCN, and CH 3 OH on 2011 September 17-21 (Rh = 2.00 
AU). Finally, on 2011 October 13 (Rh = 2.00 ALT), DiSanti et al. (2014) obtained a (poorly 
constrained) rotational temperature for water T rot = 421 2 ]] K. 


Also the column densities for CH 3 OH, have been derived from the rotational diagram 
and are listed in Table 2. 


4.2. Column densities, Production rates and Abundances ratios 

Beam-averaged column densities were derived for CO, HCN, CS and o-H 2 CO assuming 
that the cometary coma filled the beams of the respective telescopes. The column densities 
for the observations were calculated using equation 2 from Gicquel et al. (2014) with T rot = 
50 K. The column densities derived for CO, HCN, CS and o-H 2 CO are listed in Table 2. 




E„/k [K] E u /k [K] 


Fig. 2. — Detection of CH 3 OH (J=5-4) taken on 2012 February 19 (o), and March 26 (p) 
with the SMT and rotational diagram. See Table 1 for dates, parameters of observations, 
and reference letters on the right-hand side of each panel. (*) denote CH 3 OH lines. 

The production rates were derived for CO, HCN, CS and o-H 2 CO from a Monte Carlo 
model described in Milam et al. (2004, 2006). CO and HCN, are calculated as parent 
molecules. CS is considered as both a parent and a daughter species. o-H 2 CO 
is calculated as both a parent and an extended source species. In the case of CO, 
HCN, CS and o-H 2 CO as parents molecules, we used the velocity law v p = 0.8 R^°' 5 [kms -1 ] 
from Biver et al. (1999). The determination of the velocity from the CO, CS and 
o-H 2 CO lines are difficult so the integrated fine intensity is derived directly from 
a Gaussian fit (e.g. section 3). To confirm the hypothesis of the above law, 
we derived velocity FHWM from fits of the HCN lines. The velocity derived 
from the observations are in good agreement with the velocity law. To conclude, 
for simplification and coherence in the analysis we adopt v p defined above for all 
observations. For o-H 2 CO as an extended source, we used parent scale lengths in the range 
(4000-8000)7?^ 5 [km] (after Biver et al. (1999)) and a velocity law equal to v p = 0.6 i ?]) 0 ' 25 [km 
s -1 ] for the parent species (Combi & Fink 1997; Milam et al. 2006) and Vd = 0.8i ?^ 0 -5 [kms -1 ] 
for the daughters molecules. From the first observations of comets with ALMA, Cordiner et 
al. (2014) found smaller scale lengths (a few thousand kilometers in comet Lemmon (at Rh 
= 1.5 AU) and only a few hundred kilometers in comet ISON (at R^ = 0.45 AU) probing 
only the innermost few thousand km of the coma. Our observations have distances up to ~ 
20,000 km from the nucleus, so the scale lengths of Biver et al. (1999), derived from single 
dish observations, are more applicable for these data. For CS as a daughter molecule, 
we assumed a production from a short-lived species such as CS 2 (Snyder et al. 





2001) with a lifetime equal 370 Rf [s] (Huebner et al. 1992). CS 2 has much smaller 
photo-dissociation scale length 300 km) than the beam radius (1.5 x 10 4 km), 
so the difference between the gas production of CS as daughter molecule or as 
a parent molecule is negligible (DiSanti et al. 2009). 

Table 2 summarizes the production rates of the molecules observed toward comet C/2009 
PI (Garradd). This model assumes isotropic outgassing, which is reasonable for the analysis 
of these data since the model simulates the observed column densities within a large beam 
(with respect to the comet). The uncertainties introduced into the modeled production rates 
are dominated by the errors on the measured column densities. We used the production rate 
of water from Combi et al. (2013), Feaga et al. (2014) (on 2012 March 26), Bockelee-Morvan 
et al. (2014) (on 2012 February 17 and 23) and Bodewits et al. (2014) (on 2012 April 24) to 
obtain the Q/Qh 2 o ratio with the SMT and JCMT. 

Combi et al. (2013) observed water every 2 days on average from 2011 August 
15 to 2012 April 6 with SWAN (Solar Wind Anisotropies). Also, Bodewits et al. 
(2014) acquired 107 observations from 2011 April to 2012 October with Swift/ UVOT. They 
both identified an asymmetry of Qh 2 o about perihelion, with larger values before perihelion 
than after and different trends with heliocentric distance. The gas production rate of water 
peaks around 50 - 100 days before the perihelion (Figure 3). In this figure, we used the 
curve from Bodewits et al. (2014) to give an idea of the temporal production 
rate of water by averaging data in IR, mm-sub and UV. These results highlight the 
importance of using the Qh 2 o at the same date as our observations if the correct abundance 
ratio is to be determined for comparisons to other targets and dates. 



Table 2 . Beam average column densities, photodissociation rates, production rates of HCN, CO, CS, 0-H2CO and 
CH3OH, production rates of water and the ratio Q/Qh 2 o with the SMT and JCMT. 
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l The photodissociation rates are from Photo Ionization/Dissociation rates (Huebner et al. 1992)). 
* Qh^O from Combi et al. (2013) at comparable dates. 

Q H 2 O from Bockelee-Morvan et al. (2014). 



*Qh 2 0 from Feaga et al. (2014). 
Qh 2 0 from Bodewits et al. (2014). 
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From our SMT and JCMT data we obtain the following average abundance ratios, 
relative to water: CO: 7.03 ± 1.84 %, HCN: 0.04 ± 0.01 %, o-H 2 CO: 0.14 ± 0.03 % as 
a parent molecule (and 0.28 ± 0.06 % as an extended source), CS: 0.03 ± 0.01 % (as 

a parent or daughter molecule) and CH 3 OH: 3.11^0 5 ® %. The production rate of 
CS as a parent or daughter molecule are similar due to the short lifetime of the 
parent specie. 


4 . 2 . 1 . CO 

At mm-submm wavelengths with the IRAM-30m, Biver et al. (2012) found a significant 
difference in the abundance of CO relative to water between 2011 October 13-21 (5 %) and 
2012 February 15-19 (10 %). At subnnn wavelengths with the JCMT, comet Garradd was 
observed in four observing runs between 2011 July and 2012 January, and Yang & Drahus 
(2012) deduced a preliminary C0/H 2 0 = 7.02 % on 2011 September 23-25. The only CO 
abundance ratio from observations in the ultraviolet (~ 20 %), was derived by Feldman et 
al. (2012) using the Hubble Space Telescope on 2012 January 19. In the infrared Paganini 
et al. (2012), Villanueva et al. (2012), DiSanti et al. (2014) and Feaga et al. (2014) reported 
respectively 12.5 %, 10.4 ± 1.5 %, 9.12 ± 1.22 % and 60 ± 1.22 %. (Paganini et al. (2012) 
acquired spectra of CO with the CRIRES at ESO’s Very Large Telescope on 2011 September 
18 and 21; Villanueva et al. (2012) detected CO on 2011 September 8 with the CSHELL 
at IRTF; DiSanti et al. (2014) conducted observations of CO on 2011 October 13 with 
the NIRSPEC at the Keck II; Feaga et al. (2014) measured CO with the High Resolution 
Instrument Infrared Spectrometer (HRI-IR) on board the Deep Impact Flyby spacecraft on 
2012 March 26 and 2012 April 02). Our value of CO/H 2 0 (7.03 ± 1.84 %) is consistent 
with all the values determined at mm-submm and IR except the one reported by Feaga et 
al. (2014). Indeed, Feaga et al. (2014) show that Garradd exhibited a monotonic increase in 
CO production throughout the entire apparition (which seems to finally decrease with our 
observations late March (Figure 3). However, with a decrease in water after perihelion and a 
increase of CO, Feaga et al. (2014) obtained the highest CO/H 2 0 ever observed inside 3 AU 
for any comet (~60 %). Previous CO-rich comets have had CO/H 2 0 values ranging from 
10-30 %, 4 have ratios >10 %, with only one, C/2008 Q3 (Garradd), reaching « 30 % (at 
R h = 1.7 AU). On 2012 March 26, our value for the gas production rate of CO is 
lower than the one derived by Feaga et al. (2014). By using Qh 2 o from Feaga et 
al. (2014) instead of Combi et al. (2013), we derived a abundance ratio 3 times 
lower (Table 2). Feaga et al. (2014) explained the CO-rich composition in some comets 
by the possible formation of CO in the outer regions of the disk (where the stellar radiation 
is less intense) or by capture of CO within the water ice mantle. Due to the different 
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abundance ratio bu using Qh 2 o from Feaga et al. (2014) or Combi et al. (2013), we 
used also used Qh 2 o from Bockelee-Morvan et al. (2014). We conclude that the 
difference between the two abundance ratio on 2012 February 18 is negligeable 
(Table 2). 


4.2.2. HCN 

We derived abundance ratio of HCN relative to water between (3.52 ± 1.53) x 10~ 3 
% and (6.77 ± 0.57) x 10~ 2 %. By averaging all the data we obtained a surprising low 
abundance ratio close to 0.04 ± 0.01 %. Preliminary average ratios have been obtained at 
mm-submm wavelengths with the IRAM-30m and the JCMT. Both Biver et al. (2012) and 
Yang & Drahus (2012) derived an average HCN/H 2 0 ~ 0.11 % from data before and after 
the perihelion. In the infrared Paganini et al. (2012), Villanueva et al. (2012), DiSanti et al. 
(2014) found abundance ratio between 0.19 ± 0.02 (with the NIRSPEC at the Keck II on 
2012 January 8) and 0.37 ± 0.05 % (the CRIRES at ESO’s Very Large Telescope on 2011 
September). As explain in section 4.2.1, Biver et al. (2012) found a significant difference in 
the abundance ratio before and after perihelion. Considering the variation with time, the 
only comparable value after perihelion in given by DiSanti et al. (2014) on 2012 January 
8. However, our values with the SMT are lower (2011 December 29, 2012 January 18 - 19) 
or around (2011 December 28, in the error bars) the typical factor of 2 between mm- 
submm and IR measurements (Villanueva et al. 2013). As Feaga et al. (2014), we supposed 
a seasonal effect or a characteristic of a relatively young comet, but we didn’t derive the 
rotational period due to the lack of time resolution data. We obtained the minimum HCN 
ratio relative ((3.52 ± 1.53) x 10 -3 %) on 2012 April 24. However, we used the Qh 2 o at the 
closest date (2012 April 30) from Bodewits et al. (2014) due to the lack of complimentary 
observations on 2012 April 24. The low ratio of HCN relative to water can be explained 
by uncertainties in the water production rate, and possibly the FOV of observations made 
with different techniques. To investigate it, we derived the abundance ratio on 2012 
March 26 from Combi et al. (2013) (UV) and Feaga et al. (2014) (IR). We found 
a abundance ratio 2 times higher with the Qh 2 o from the IR data (Table 2). 
However, the difference between the abundance ratio with the Qh 2 o from Combi 
et al. (2013) or Bockelee-Morvan et al. (2014) on 2012 February 18 and 23 is 
negligeable (Table 2). 
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12.3. CS 

Our derived abundance ratio of CS relative to water (0.03 ± 0.01 %) is in good agree- 
ment with the average values derived by Biver et al. (2012) using the IRAM-30m (0.04 %), 
but by Yang & Drahus (2012) in 2011 September with the JCMT (0.07 %). As explained 
in section 4.2, the production rate of CS from the nucleus or from a short-lived 
species such as CS 2 are similar. We derived mixing ratio on 2012 February 18 
and 23 by using Qh 2 o from Combi et al. (2013) and Bockelee-Morvan et al. (2014) 
and found similar results (Table 2). 


4 . 2 . 4 . h 2 co 

The relative abundance of o-H 2 CO as an extended source (0.28 ± 0.06 %) is consistent 
with the average value from Biver et al. (2012), (0.30 %, from submm observations), while its 
abundance as a parent molecule (0.14 ± 0.03 %) is consistent with that derived by DiSanti 
et al. (2014) (0.11 ± 0.041 %), obtained on 2011 October 13 with the NIRSPEC at the Keck 
11. The held of view in the infrared is typically just a few arcseconds, so extended H 2 CO 
sources may not be as obvious in those observations. 


4-2.5. Methanol 

Our estimate for the abundance ratio of the methanol relative to water, 3. Ill,)'®® %, 
agrees with those by Bockelee-Morvan et al. (2014), (3.40 ± 0.6%, using the Herschel Space 
Observatory on 2012 February 22) and DiSanti et al. (2014), (2.14 ± 0.34 %), although all 
these measures are notably larger than those by Biver et al. (2012) (1.6 %), Yang & Drahus 
(2012) (1.85 %), and lower than the one by Paganini et al. (2012), (3.90 %). 


4.3. Temporal evolution 

Our resulting production rates of HCN and CO and others from the literature are plotted 
with time and heliocentric distance in Figure 3. Observations of comet Garradd spanned 
4 months for CO and HCN post perihelion. The temporal evolution of CO and HCN are 
available in the online journal as supplemental information. 

As explained in section 4.2, Combi et al. (2013) and Bodewits et al. (2014) studied in 
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detail the temporal evolution of water. The water production rate increased as a function of 
i?/ 6 then peaked around 50-100 days before the perihelion (Feaga et al. 2014; Bodewits et 
al. 2014). Qh 2 o remained constant for about 100 days, then decreased as a function of R^ 4 . 
Both Combi et ah (2013) and Bodewits et ah (2014) also showed that estimates of the water 
production rate increased with the observed field of view (FOV). Combi et ah (2013) and 
Bodewits et ah (2014) thus supported the suggestion of Paganini et ah (2012), Villanueva 
et ah (2012) and Bockelee-Morvan et ah (2012) that a significant fraction of the water was 
from an extended source of icy grains, especially pre-perihelion. As suggested by DiSanti et 
ah (2014) the difference between pre- and post-perihelion might also be due to the seasonal 
variation in the heating of one or more active region of the nucleus. 

In the case of CO, the gas production rate after perihelion remains essentially constant - 
or perhaps actually increases slightly: the CO production reported by Feaga et ah (2014), 97 
days after perihelion, was three times higher than our value obtained 94 days after perihelion. 
Such behavior for CO, constant or monotonically increasing, has never been observed before. 
Feaga et ah (2014) suggest that such striking asymmetric activity is either a seasonal effect 
or a characteristic of a relatively young comet. This could also be a minor outburst event of 
highly volatile species. 

In the case of HCN, the gas production rate increases before perihelion - or perhaps 
remains constant and decreases after perihelion. The temporal evolution of HCN tends to 
be similar to the one of water. However, the lack of data does not allow us to study the 
behavior the gas production between 70 days before perihelion and the perihelion. Also, 
the temporal evolution of the gas production rate, of the line profiles or of the 
shift of the peak line with the SMT shows variation that may be attributed to 
the rotation of the nucleus. We tried to obtain the rotational period of the nucleus but 
the lack of time resolved data made it impossible. 


5. Conclusions 

We conducted observations of comet C/2009 PI (Garradd) from 2011 December 28 to 
2012 April 04 using the SMT and the JCMT and have reported here detections of HCN, o- 
H 2 CO, CS, CO and CH 3 OH. We have derived column densities, production rates, rotational 
temperature and relative abundances. We concluded that Garradd is normal in CH3OH, 
depleted in HCN, o-H 2 CO and CS and slightly enriched in CO, which is in good agreement 
with other studies. Temporal variations in Garradd’s HCN are similar to water (but with 
short-term variation), with a decrease after perihelion, while that of CO in the same period 
shows an unexpected constancy or even increase. The observed short-term variability of HCN 
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Fig. 3. — Comparison of the evolution of the production rate for HCN and CO around comet 
Garradd’s 2011 perihelion with the H 2 0 curve. Bodewits et al. (2014) average data 
from Combi et al. (2013), Paganini et al. (2012), Bockelee-Morvan et al. (2012), 
Bockelee-Morvan et al. (2014), DiSanti et al. (2014), Villanueva et al. (2012), 
Feaga et al. (2014) and Feldman et al. (2012). 

may be due to a seasonal variation or to Garradd’s relative youth. These results highlight 
the importance of studying young comets to better understand their temporal evolution. 
Long term evolution, complimentary studies, and time resolved data are important for the 
full analysis and interpretation of cometary data. 
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